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We have found a new class of time dependent partial waves which are solutions of 
time dependent Schrodinger equation for three dimensional harmonic oscillator. We 
also showed the decomposition of coherent states of harmonic oscillator into these 
partial waves. This decomposition appears particularly convenient for a description of 
the dynamics of a wave packet representing a particle with spin when the spin-orbit 
interaction is present in the hamiltonian. An example of an evolution of a localized 
wave packet into a torus and backwards, for a particular initial conditions is analysed 
in analytical terms and shown with a computer graphics. 



PACS number(s): 03.65.Ge, 32.90+a 

1 Introduction 

The rapid technological development of short-pulsed lasers during the last decade made 
it possible to produce and detect particular states, coherent superpositions of stationary 
electron states for a wide variety of physical systems. The dynamics of these initially well 
localized wave packets is a subject of much current investigation in many areas of physics 
and chemistry |l|, [|. An extensive research by means of both theoretical and experimen- 
tal methods has brought an understanding of the intriguing phenomena of a hierarchy of 
collapses and revivals for particularly prepared wave packets in Jaynes-Cummings Model 
(JCM) of quantum optics [3-7] and in Rydberg atoms [8-14]. 

In present paper we construct new solutions of time dependent Schrodinger equation for 
spherical harmonic oscillator (HO). We use these solutions further for description of the 
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evolution of wave packets representing a particle with spin moving in HO potential with 
additional spin-orbit interaction. Time dependent partial waves allow for much clearer in- 
terpretation of subtle interference effects as well as for a substantial acceleration of numerical 
codes used for graphical presentation of a complicated wave packet motion. 

The paper is organized as follows. In section 2 we construct time dependent partial 
waves for the harmonic oscillator and show the decomposition of coherent states into these 
states. In section 3 we show the motion of individual partial waves. In section 4 we discuss 
time dependent spinors corresponding to solutions of time dependent Schrodinger equation 
for hamiltonian containing spin-orbit interaction and show a possibility of evolution of an 
initially localized wave packet into a toroidal shape and backwards. Sections 5 and 6 present 
graphical illustrations of the wave packet motion and conclusions. 

2 Time dependent partial waves for the harmonic os- 
cillator 

It is not so well known that there exist simple time dependent partial waves which are 
solutions of the Schrodinger equation of the three dimensional oscillator. The proof will be 
given below. 

We look for solutions il)™(r, t) of the time dependent Schrodinger equation for the spher- 
ical harmonic oscillator (with h = m = uo = 1) with the following form 

^ m (f, t) = F(t) e-2 r ' 2 Wi(R e- lt r) Y™(0, if) (1) 

where F(t) and Wi must be determined in terms of a complex number R . 
We find easily that W\ must obey the differential equation 

9 d 2 Wi dWi ... . , % dF. 9n 

with a new variable Z 

Z = R e- u r . (3) 

Wi depends only on this variable if F solves the equation 

3 _ *£ = *..-. (4) 

Within an arbitrary constant factor, which is not written the solution is 

F(*)=e-§ it e-* B ° ae " a ". (5) 

The equation solved by Wi becomes 



which is the equation for the modified spherical Bessel functions. 



(6) 



In the following we will use only the modified spherical Bessel functions of the first kind 



with the usual conventions of the literature 15 
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The interpretation of these waves is very simple: there is an harmonic motion in the radial 
part of the wave function keeping the angular part the same because of angular momentum 
conservation. 

The partial waves ^ m (r, £) just defined occur in a natural way in the expansion of a co- 
herent time dependent gaussian wave packet into partial waves. In such a wave packet the 
constant Rq finds its interpretation by combining the position of the center of the wave 
packet r*o to its mean momentum p Q at a time t = 0. Let r t and p t be such vectors at time t 
and let us define 

Rt = f t + ip t , (9) 



while 



These two vectors are related by 



Ro=r + ip . 



(10) 



Rt = R e- u . (11) 

A normalized gaussian wave packet centered on r* with mean momentum p t is now written 



as 



^(r) = vr-t e -W-*tf e ipt * 



-2 -I r 2 -±rt 2 r-Rt 

ix 4 e 2' e 2 ' e * 



(12) 



Using 10.2.36 from reference |15[ the modified Bessel function of the first kind of argument 
Z = R t r = Rq e~ lt r appears in the expansion of the last exponential of (12) 



-f-Rt 
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and Rn is defined as 



<■■- ="£(21 + 1) J— WZ)P z (cos9) 



Ro = (Rq ■ Rq)2 = (r 2 - p 2 + 2ir ■ p ) 2 



(13) 



(14) 



In a gaussian wave packet all the time dependent partial waves share the same parameter Rq. 
However it is possible to consider more general wave packet for which Rq might be different 
for different I. 

The Legendre polynomial -P^cos 0) depends on the (complex) angle G between f and R t . 
However this angle is time independent. Indeed for each direction a, a = x,y,z one has 



r a (t)+ip a (t) = (r«(0) +i Pa (0))e- u 
as well as equation (11), therefore 

x[x{t) + ip x {t)] + y[y(t) + ipy{t)} + z[z(t) + ip z {t)] 



cos(9) 



Rtr 

x[x(0) + ip x (Q)} + y[y(0) + ip y (0)] + z[z(0) + i Pi 

R r 



(15) 



(16) 



Then if the complex direction 9 Ro , (p Ro of R is introduced through 

z(0)+ip z (0) 



cos 9 



Ro 
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one writes the addition formula: 



. x{0)+ip x (0) 

sm9 Ro cos^p Ro = (17) 

. . ■ y(o) + ip y (o) 

sm9 Ro smtp Ro = - 
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Pi(cosG) = — — £ (-irYr(9^)Y l - m (9 Ro ^ Ro ) . (18) 

rn=—l 

We can now rewrite ip^ as 

oo I 

$sjf) = 4^e-i(^?) £ E (-l) m ^r m (^oW^^)in^) . (19) 

Z=0 m=-l 

We will now extract the function F (5) in order to show explicitly the time dependent partial 
waves 



— e 2 r * 6 2 rt 'P l g 2 

The classical energy i? = r ° 2 ? ' has been introduced above. The time dependent partial 
waves (1) which appear in the gaussian wave packet at time t are then expressed as: 

2 2 

?/f (f,£) = e -('iit+^fPt+ r -^-^) w l (R t r)Y l m (9,ip) . (21) 

The gaussian wave packet (19) can now be expressed as follows in terms of (21) 

i> A (r) = An^ e^ t+r ^-^Y,(-^) m Yr m (^^W(r,t) . (22) 



I irS. f i n-p t \ go 
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Clearly one obtains a solution of the Schrodinger equation only if we change the phase. Let 
ip Rt (r) be this solution 

^R t ( f ) = e^ ( t t+ ^ } ^ t (f) 

= Y,(- i r^^ Y r m (°Ro^R )^r(r,t) (23) 

tin 

lm 

The weight Ci m (Ro) of the time dependent partial wave tpTi^-i t) * s now seen to be 

Q m (Ro) = (-l) m ^e-^Yr m (9 RoJ ip Ro ) . (24) 



An expansion similar to equation (23) is found in reference \W\. However the time 
dependence of the radial function is lacking in this reference as well as the phase which is 
needed to correct the gaussian wave packet. The existence of this phase is mentioned in 



some textbooks [17 . 



3 Representation of the time dependent partial waves 

Let us define the phase of i?o by 

Ro = (ro 2 - Po 2 + m •#>)* = \Ro\ e lSo . (25) 

The argument of the modified Bessel function is Z = r|i?o|e _J( ' t ~ <5o * ) . Changing the initial 
conditions of the wave packet produces a change of |i?o| and of 5q. However it is clear 
that the change of the phase can be taken into account in the radial partial wave Wi(z) by 
a shift of the origin of time. It is therefore possible to analyse the different partial waves 
by assuming that So = 0. This possibility arises if po = 0, Ro = r Q . For these values the 
gaussian wave packet performs a linear motion. It is striking that the radial motions of the 
partial waves corresponding to such a case contain all the proper information which can be 
used for all the different trajectories, i.e. circular or elliptical. The differences between the 
trajectories will be introduced by the coefficients Ci m (Ro). 

In the following figures we have chosen to study the radial waves for the mean energy 



E Q = N = 20 i.e. r = \ / 2E = V40 = R . The 8 lowest partial waves present in the 
development of a gaussian wave packet are represented in figure 1 [as a matter of fact we 
have represented the density multiplied by r 2 ] for a time range going from t — to |- (T 
is the harmonic oscillator period). The densities exhibit the same motion toward the origin 
and of course this radial motion is symmetric with respect to t — j. However there is clearly 
an effect of the centrifugal barrier since the waves with the highest / are repelled from the 
origin. Also some secondary maxima are present. It is worth to mention once again that the 
effects of the centrifugal barrier cancel when one adds the waves to produce a linear motion. 
We should also note that the spread of the partial waves does not depend on / in a significant 
manner. 

In figure 2 the squares of 16 lowest partial waves (\ripi(r, t)\ 2 ) are represented as functions 
of r for t — 0. The concentration of the waves at the same r is spectacular. The intensity 
of the waves increases when I goes from I = to I = 4 and decrease afterwards. In order to 
explain the concentration of the wave in three dimensions the angular part plays a role as 
will be discussed later on. 



4 Time dependent spinors 

We will now study the evolution of a coherent wave packet which is initially in a pure spin 
eigenstate chosen as s z = +|. This wave packet will evolve on the action of the harmonic 
oscillator evolution operator Uo(t) and of the spin-orbit part Ui s (t). As in our previous 
papers [18-22] we will use the spin-orbit hamiltonian with a constant factor k: 

V 1s = k (f • a) (26) 

and therefore U$ commute with VJ S . The hamiltonian contains the ratio of the two time 
scales (let us recall that uj — 1) 

Ti s 2ix uj 1 ,„„. 

-£ = — 7T = - ■ ( 27 ) 

1 K Z7T K 



The fully normalized spinor is written at time t = as 

i s jr^[\^)^c lmi s. 0) (^my (28) 

We have shown previously [4] that 

U ls (t) = f(t)+g(t)(l-a) (29) 

with the following operators 

f(t) = e *l(co S ^-± s infi~) (30) 

g(t) = ^(-|sinfi0 (31) 



n\lm) = yi + 4l 2 \lm) = (2l + l)\lm) . (32) 

Uo(t) transforms Rq into R t while Ui s (t) acts on the spherical harmonics present in the 
ip™ (r , t) . The result is 



^ t (r) = e H(§^)-^-^] (33) 



where fi(t) and gi(t) are simply obtained from (|30| ) and ( pT|) by replacing $1 by (2/ + 1). This 
expression corresponds to the most general initial condition with the assumed initial spin 
direction. The presence of the time dependence in the functions / and g from one hand, the 
presence of Y" ; m+1 in the lower part of the spinor are the two differences which change in an 
appreciable manner the time evolution in the case where k^O. 

We want now to concentrate our efforts on the evolution of a wave packet which is 
cylindrically symmetric around Oz for t = 0. We assume therefore (denoting unit vector in 
Oz direction by z) 

r = -r z , p = , Ro = r , (34) 

2 

£o = y , 0^ = n , <^ = , (35) 

C lm (R ) = 5 m0 (-1) 1 2vr-3 e -f y/2l + 1 . (36) 

Let us study the sign of the product Czo(-Ro) Y®{9, ip) along the Oz axis i.e. 9 = it and 9 = 0. 
One has 



y,>,°) = (-i)'v^ < 37 > 



y <>- 0) = V^' (38) 



All the partial waves act coherently in the direction 9 = n while there is a destructive 
interference for 9 = 0. We also know that there is a radial concentration of the wave packets 
shown on figure 2. This explains mainly the known result that the wave packet is initially 
a gaussian centered at Fq = — To z. 



In the direction 9 = | we have 



Y l °(^,0) = if /is odd (39) 



; 



while for / even its sign is (—1)2. A destructive interference between the even states is also 
predicted on the xOy plane with a different character as for the +Oz direction. 

We will now try to repeat these arguments for t ^ for each component of the spinor 
(33). For the part with s z + \ we must find the phase of the product CiqY^O, <p) fi(t). At 



a time t — -^ = | (assuming k — 1) we have for high I: 

/,(£) » e'* cos(2/ + if- = e *5 ^ (-l)i . (40) 

The coherence of the product considered is now obtained for 9 — | and for even Z. It is 
interesting to see that the part with s z = +| has still cylindrical symmetry around Oz but 
moreover if t = -}£■ or we shall have the result of figure 2 i.e. a high radial concentration 
of the wave. One readily understand that at this time the wave is highly peaked on a vortex 
of radius tq with symmetry around Oz, the smaller radius of the vortex being roughly the 
initial radial spread. 

As for the part with s z = —\ we must discuss the sign of the product C/o Y l 1 (9, p) giif) 
x Jl{l + 1) also for t — | and for 9 — | . Now 1^ 1 (|, (/?) = if / is even, its sign for odd I is 
(— 1)~5~ while <7j(|) has also this sign since 

.7r, ,t 2i . , . n —2i , . i-i -7t v2 

o,(-) = -e*4 _ S in 2 + 1 - = (-1)~ e l * — . (41) 

y v 2 ; 2/ + 1 v ! A 21 + 1 V ; 2 v ; 

However we have also Y^(9,(p) rs e Jip . In the case of the part with s z = — | it is the odd 
partial waves which produce a vortex similar to that to the part with s z = +| in the limit 
of high /. A difference occurs in a tp dependent phase. 

For intermediate values of the time it is not possible to provide a similar discussion. 
However, by continuity we can understand that the vortex rings are created from the spherical 
initial wave packet and get their full extension for the configuration and at the time chosen. 



5 Numerical simulation of the vortex rings 

The manifestation of the vortex rings depends on the parameter k. It is simpler to freeze 
the evolution under Uq and to consider only the spin-orbit evolution to begin with in order 
to emphasize the effect. Figures 3 and 4 represent the density of the full wave packet at 
various times on the planes yOz and xOz, respectively. Here also we have considered the 
wave packet which was analyzed in figure 2. [For convenience the spin direction is s x = +| 



and r = r x at time t — 0]. At time t = a gaussian is represented. At time £ = -^ a 
vortex ring has been created and gets the maximum radius at t = ^ when intersecting the 
plane yOz. Moreover figure 4 shows that even a second vortex with a smaller amplitudes 
has also been created that we cannot explain into simple terms. At time t = ^^- the vortex 
ring has a decreasing radius and is centered on a point with x < 0. Finally at time t — ^ 
the wave packet is reassembling near a point with r = —r Q x. At this time we have shown 
in our previous papers that the spin is approximately reversed. Finally the behaviour for 
Ti s /2 < t < Tis has not been shown for it is totally reversible if we assume a frozen oscillator. 

The vortex rings evolve approximately on a sphere. The intersections of the wave packet 
with a plane xOz shown on figure 5 present indeed the feature that the distance to the center 
is about ro at all times. 

The following figure 6 shows the decomposition of the toroidal wave packet into its spin 
components in a direction perpendicular to the classical motion. The left and central columns 
display the shapes of subpackets in which spin field is antiparallel and parallel to em Ox axis, 
respectively. Due to symmetry with respect to em Oz axis (classical trajectory), the similar 
decomposition with respect to spin components in an arbitrary direction perpendicular to 
em Oz axis results in a picture rotated by a certain angle. 

6 Conclusions 

The dynamics of wave packets in static potentials implies subtle interference effects which 
have been beautifully illustrated in the case of Coulomb potential [9-14,16]. Two general 
mechanisms have been identified: one is the spread of the wave on the top of a classical 
trajectory, the other is the regime of partial and also almost complete recurrences. In the case 
of a pure harmonic oscillator the existence of a single frequency allows a unique interference 
mechanism of the partial wave that leads to a dispersionless coherent wave packet. This is 
the case where quantum mechanics comes closer to classical mechanics: the dynamics of the 
density probability is identical to the dynamics of the density distribution of the classical 
ensemble. This property is expected to be lost if a perturbation is added to the potential. If 
the perturbation is a spin-orbit potential the dynamics of the problem resembles very much 
that of the JCM |24]] . There exist two versions |25[ of this model where the time evolution is 



exactly periodic in the closest analogy to our model: the Raman coupled model and the two 

photons JCM. In the harmonic oscillator with spin-orbit the dynamics is however richer if 

one follows the angular evolution with different initial conditions. We have shown previously 

T9fl that the amount of partial revival for half a spin-orbit period depends indeed on these 



conditions. We have discussed extensively [|18|, |21| cases where the wave packet is divided 
into two parts which rotate in opposite directions on a circular average 

trajectory in much the same way as the Stern-Gerlach effect. In the present paper we 
have shown a new dynamical behaviour for cylindrically symmetric initial conditions: vortex 
rings are created and destroyed periodically. In a way this is a case where the coherence of 
the wave is kept for the r and 9 coordinates but where the spread is applied only to the 
variable. It is an open question whether this effect relies only on the properties of the partial 
waves of the harmonic oscillator discussed in this paper. 
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Figure captions 

Fig. 1. Radial motion of eight lowest partial partial waves with I — 0, 1, . . . , 7 contributing 
to the linear motion of the gaussian wave packet. Shown is |r^(r, t)\ 2 as a function of 



T 

vertical scale is the same for all figures. 



r for 30 time steps in the interval t G [0, -^]. The case Eq = N = 20 is presented. The 



Fig. 2. Decomposition of the gaussian wave packet into partial waves at t — 0, tq — y2N = 
\/40. The lowest 16 partial waves with / = 0, 1, . . . , 15 are shown. Waves with / = — 4 
are represented by solid lines (intensity increases with I), those with I = 5—15 by dashed 
ones (their intensities decrease). 

Fig. 3. Time evolution of the wave packet with spin under Ui s (t) operator only (evolution 
according to Uq is frozen). Shown is the |\l/(t)| 2 = |\l/ + (t)| 2 + [^.(t)) 2 as a function of 
coordinates on the plane xOz. Case N = 20. Note different vertical scales. Cases em 
a,b,c,d,e,f correspond to tj = 0, 1/8, 2/8, 3/8 15/32 and 4/8 T} s respectively. 

Fig. 4. The same wave packet as in previous figure. Here shown are cuts through planes 
perpendicular to the classical trajectory (em Oz axis ) with Zi = zq costj showing 
explicitly a toroidal shape. Time instants as in previous figure. 

Fig. 5. Contour plots of the same evolution as in figure 4, showing that centers of the vortex 
rings evolve approximately on a sphere with the radius r$. Here packets corresponding 
to different time instants are collected in the same picture. 

Fig. 6. Contour plots showing the decomposition of the total wave packet (right column) 
into parts with the opposite spin field (left and central columns). Case N — A. Top 



row a) corresponds to t — 0, bottom e) to t — \Ti 3 . Time steps are At = \Ti 
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Figure 1. Radial motion of eight lowest partial partial waves with / = 0, 1, . . .,7 
contributing to the linear motion of the gaussian wave packet. Shown is \ripi(r, t)\ 2 as 
a function of r for 30 time steps in the interval t £ [0, ?■]• The case Eq = N = 20 is 
presented. The vertical scale is the same for all figures. 
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Figure 2. Decomposition of the gaussian wave packet into partial waves at t = 0, ro = 
\/ c 2N = V40- The lowest 16 partial waves with / = 0,1,...,15 are shown. Waves with 
/ = — 4 are represented by solid lines (intensity increases with /), those with / = 5 — 15 
by dashed ones (their intensities decrease). 
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Figure 3. Time evolution of the wave packet with spin under Ui s (t) operator only 
(evolution according to Uo is frozen). Shown is the ^(t)! 2 = |\?+(t)| 2 + ^-(t)! 2 as a 
function of coordinates on the plane xOz. Case N = 20. Note different vertical scales. 
Cases a, b, c,d,e,f correspond to t 8 - = 0, 1/8, 2/8, 3/8 15/32 and 4/8 T; s respectively. 
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Figure 4. The same wave packet as in previous figure. Here shown are cuts through 
planes perpendicular to the classical trajectory (Oz axis ) with Z{ = zq cos ti showing 
explicitly a toroidal shape. Time instants as in previous figure. 
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Figure 5. Contour plots of the same evolution as in figure 4, showing that centers 
of the vortex rings evolve approximately on a sphere with the radius 7*0. Here packets 
corresponding to different time instants are collected in the same picture. 
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Figure 6. Contour plots showing the decomposition of the total wave packet (right 
column) into parts with the opposite spin field (left and central columns). Case N = 4. 
Top row a) corresponds to t = 0, bottom e) to t = i2] s . Time steps are At = ^7] s . 



